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TECHNICAL NOTE 2916 

EFFECT OF THERMAL PROPERTIES ON 
LAMINAR-BOUNDARY-LAYER CHARACTERISTICS 
By E. B. KZLunker and F. Edward McLean 

SUMMARY 



An iteration method is presented for solving the laminar -houndary- 
layer equations for compressible flow in the absence of a pressure 
gradient wherein the temperature variation of all the fluid thermal 
properties is considered. Friction and heat-transfer characteristics 
have been calculated for a stream temperature of - 67 ° F for Mach numbers 
from 1 to 10 with the use of values of the heat capacity, conductivity, 
and viscosity determined from experiment . Consideration of the temper- 
ature variation of all the fluid thermal properties causes the recovery 
factor to decrease substantially with increasing Mach number. Moreover, 
the heat -transfer rate is found to be proportional to the difference 
between an effective enthalpy, which is a function of both the surface 
temperature and stream Mach number, and the surface enthalpy. In con- 
trast, the heat-transfer rate is approximately proportional to the dif- 
ference between the recovery enthalpy and the surface enthalpy for 
solutions which employ a constant Prandtl number. The calculated skin 
friction and heat-transfer rates based upon the use of the Sutherland 
equation for viscosity and a Prandtl number of 0.75, however, are in 
excellent agreement with the results of the present analysis . 


INTRODUCTION 


Most compressible-flow boundary- layer analyses employ approximate 
relations for the temperature variation of the thermal properties . A 
simple analytic relation is usually specified for the’ variation of vis- 
cosity with temperature . Moreover, the heat capacity and the Prandtl 
number are generally considered to be independent of temperature and 
as a consequence the conductivity and viscosity are required to have 
the same temperature variation. Experiment shows that these conditions 
are not satisfied closely over the wide temperature range which occurs 
in a boundary layer in a high-speed flow. The effect of these approxi- 
mations on the calculated boundary-layer characteristics is most easily- 
studied by considering the laminar -boundary-layer flow on a flat plate . 
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Calculations of "boundary-layer characteristics which employ approxi- 
mate relations for the thermal properties are generally satisfactory as 
long as the temperature variation in the "boundary layer is not excessive; 
the range for which the solutions are satisfactory depends upon the value 
of the Prandtl n umb er used in the analysis and upon the viscosity- 
temperature relation employed. Various viscosity-temperature relations 
and values of the Prandtl number have "been used in the analysis of the 
laminar bo undar y layer on a flat plate (refs. 1 to 5, for example) . The 
choice of a Prandtl n umb er of unity (ref. 1, for example) or of a linear 
viscosity-temperature relation (ref. 2 , for example) leads to an essen- 
tial simplification, since in each case only one differential equation 
need be solved. Better solutions are obtained by taking the Prandtl 
number as approximately 0-73 (the value for air at normal temperatures) 
and by utilizing a more realistic viscosity-temperature relationship. 

The viscosity-temperature relationship is frequently taken as a simple 
power law (refs. 1 and 3, for example) or as given by the Sutherland 
equation (refs. 4 and 5, for example) . The Sutherland equation for vis- 
cosity shows especi ally good agreement with experiment over a wide temper- 
ature range. 

At high supersonic Mach numbers, where the temperature variation 
within the boundary layer is large, the assumption of a constant Prandtl 
number and the vise of an approximate viscosity-temperature relation is 
not jusitifed a priori. A few calculations of the laminar -boundary- 
layer characteristics on a flat plate have been made which have not 
imposed these restrictions . Experimental values of the thermal properties 
were utilized in some calculations of the laminar -boundary- layer charac- 
teristics by an iteration method at Mach numbers from 1 to 10 in refer- 
ence 6 and the thermal properties employed in the differential-analyzer 
solutions of references 7 and 8 were taken in part from experiment. For 
the latter two analyses, the stream conditions were determi n ed from the 
solution. Consequently, a graphical or trial-and-error procedure was 
required in the course of the solution to obtain values corresponding 
approximately to the desired stream conditions. In the iteration methods 
used in reference 6 and in the present analysis, the solutions are deter- 
mined for specified stream conditions . 

The present investigation is a continuation of the work reported in 
reference 6 . The purpose is to determine accurately some characteristics 
of the boundary layer at high supersonic speeds and to show the effect 
of various approximations to the thermal data of the fluid on the accuracy 
of the calculated friction and heat-transfer characteristics. Tabulated 
values of the thermal properties, found from experiment, have been 
employed in the calculations in order to obtain reliable data. Friction 
and heat-transfer characteristics, calculated for a stream temperature 
of - 67 ° F for Mach numbers from 1 to 10, are presented in the form of 
graphs . Tabulated values of the corresponding velocity and temperature 
profiles and their derivatives are available upon request from the 
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National Advisory Committee for Aeronautics . The iteration method 
described herein for the solution of the equations is similar to that 
presented in reference 6. The computations are facilitated however 
through a different formulation of the equations; the convergence of 
the iteration procedure is improved considerably and a troublesome 
singularity is avoided. 


SYMBOLS 


x >y 

u,v 

P 

P 

i 

k 

T 

C P 



v 

U 1 

Pr 






coordinates para ! l ei and normal to stream direction 
components of velocity along x- and y-axes 
density 

coefficient of viscosity 

enthalpy per unit mass 

thermal conductivity 

absolute temperature 

heat capacity at constant pressure 


kinematic viscosity 
stream velocity 

Prandtl number 

i - ii 

enthalpy function, - — 

u i 2 / 2 

integrating factors 



3 dTj 
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r shear stress at surface 

R Reynolds number 

Cf average skin-friction coefficient 

q local heat-transfer rate per unit area at surface 

7 ratio of heat capacities at constant pressure and con- 

stant volume 

M Mach number 

r temperature recovery factor 

c constant 

Subscripts : 


1 

0 

e 

r 

s 

Superscript : 
* 


free stream or at edge of boundary layer 

stagnation condition of stream 

effective 

recovery 

surface 

dimensionless quantity based on stream conditions 


ANALYSIS 


Basic equations .- With the x-axis taken in the free-stream direc- 
tion, the •two-dimensional compressible-flow boundary-layer equations 
for steady motion in the absence of a pressure gradient may be written 
as follows: 

Momentum equation: 


pu 


du 

Sx 


+ 


du d / chi\ 

PV Sy = M by) 


( 1 ) 
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Continuity equation: 


£ (pu) 


+ 


0 


( 2 ) 


Energy equation: 



+ pv 


3i 

dy 




(3) 


where u and v are the velocity components in the direction of the 
x- fmrl y-axes, respectively, p is the density, T is the absolute 
temperature, i -is the enthalpy per unit mass, H is the coefficient 
of viscosity, and k is the thermal conductivity. The heat-conduction 
term of equation (3) may be expressed in terms of the enthalpy. With 

the use of the relation c_ = — , where o_ is the heat capacity at 

* dT * 

constant pressure, the energy equation becomes 


pu + pv 
ox ay 


JL/JL £*\ 

^y\c p by) 




For a flow with heat transfer (specified constant surface temper- 
ature) or for a flow with zero heat transfer, the three partial- 
differential equations which describe the motion in the boundary layer 
may be reduced to a set of two ordinary differential equations . For 



Up is the stream velocity and Vp is the kinematic viscosity based 

on the stream temperature , is taken as the independent variable. Trans- 
forming equation (2) and solving for pv results in 



3 -y-(pu) dq 
dq 
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and after an integration by parts 



Transforming equations (l) and (k) and making use of equation (5) gives 



With the introduction of the Prandtl number Pr, the nondimensional 
function 





and the d i mensionless quantities 


* 


u 






equations (6) and ( 7 ) become 



(8) 
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and 


dAL ii'l + Pr,f ai/id, ^ = -SPr.^/Ulf 
dq\Pr dr)/ h* \Pr dr]/ ^ \dr) / 


(9) 


where 


f (tj) = - f P*n* dr] 

2 Jo 

The subscript 1 refers to values at the edge of the boundary layer 
or values in the free stream and iQ is the stagnation enthalpy of 

the stream. 

The transformation of the boundary-layer equations employed herein 
follows essentially that of Schuh (ref. 9) * This formulation of the 
boundary-layer equations for the flat plate has certain advantages over 
those in which the velocity (refs. 4 and 5, for example) or a function 
related to the stream function (ref. 6, for example) is used as the 
independent variable . A singularity arises in the solution of the 
equations in the latter two forms which complicates the numerical solu- 
tion of the equations . In contrast, no singularity occurs in the solu- 
tion of equations (8) and ( 9 ) • Moreover, the iteration method used to 
solve these equations converges satisfactorily. 

Equations (8) and ( 9 ) can be solved numerically by a method of 
successive approximations. The thermal quantities are functions only 
of the temperature and consequently functions of the independent vari- 
able T] . Therefore, the coefficients of the derivatives in equations (8) 
and (9) may be considered as known functions of T) if some initial 
solution is given. With the coefficients considered as known functions 
of the independent variable and with the integrating factors 



a = e 

and 


0 = e 


'^f 


f * 

— Pr dT) 

0 H 
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the equation of motion (8) and the energy equation (9) may "be written as 



Solution of boundary- value problem .- The solution of equations (10) 
and (ll) gives the distribution of velocity and enthalpy throughout the 
"boundary layer subject to certain conditions on the "boundaries . Solu- 
tions which have zero heat transfer at the surface and solutions with 
heat transfer (specified constant surface temperature) are considered 
herein. The boundary conditions to be satisfied are: 

For an insulated surface: 

At the surface (tj = 0) , u* = 0 and 


'8t\ 


^y / 




= Constantf — ) = 0 


y=0 


W 


T )=0 


At infinity (tj = ») , u* = 1 and -3 =0 
For a specified constant surface temperature: 

At the surface (tj = 0), u* = 0 and 3 = 3 g = Constant 
At infinity (tj = 00 ) f u* = 1 and 3=0 
The subscript s denotes the values at the surface . 

The first integral of equation (10) is 
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and a second integration gives ‘the velocity distribution 

u* = C lt T(Ti) 


where 




Since u* = 1 for rj = co, the function 
of the shear parameter is 



and the distribution 


p* = frU) 

d.T| j(oo) 


( 12 ) 


and the velocity distribution is 


u* = £(nl 
J(“) 


(13) 


Since the boundary conditions on the velocity are the same for both 
the boundary conditions of specified constant surface temperature and 
zero heat transfer at the surface, equations (12) and (13) are valid for 
both surface conditions. 

The temperature or enthalpy distribution throughout the boundary 
layer is obtained by integrating equation (ll) . The first integral of 
equation (ll) satisfying the condition of zero heat transfer is 


p* do 
Pr* 


= -PUMn) 


(14) 



10 


MCA TN 2916 


where 


qp(n) 



at] 


and after a second integration the enthalpy function ■d is found as 


■9 = 




0<P &T1 


(15) 


At the surface a becomes the enthalpy recovery factor 

f “ Pr* 

*r = / — W d n 

Jo M- 


The first integral of equation (ll) for flow with heat transfer 
(specified constant surface temperature) is 



— = C p - pep 
dtj d 


where 



and the f unction Cg is dependent on the surface temperature, stream 

temperature, and Mach number which is to be determined from the bou ndar y 
conditions. A second integration gives 


* = + C 2 V " U 
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where 



From the boundary condition at infinity, C 2 is determined as 


C 


2 


■d 

e 



V(») 


where 


= U(») 


With this value of C 2 


M-* d£ 
Pr* dii 


^e - ** 

V(«) 


P - P<P 


(16) 


and the enthalpy f met ion is 


- O(l) + ^-(*e - *s) U7) 

Equations ( 13 ) and (15) give the distribution of velocity and 
enthalpy (and consequently the temperature) throughout the bo undar y 
layer for the insulated surface whereas equations ( 13 ) and ( 17 ) give 
the distribution of velocity and temperature corresponding to a specified 
constant surface temperature . The distributions of shear and heat trans- 
fer are determined from equations ( 12 ) and (l 6 ) . 
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Method of numerical solution .- The solution of the equations for 
the velocity or temperature distributions is one of successive approxi- 
mations; consequently, some approximate solution for the distribution 
of velocity and temperature is required. The Blasius solution may be 
taken for the distribution of velocity and the parabolic relation between 
temperature and velocity for a Prandtl number of unity (ref. 10) may be 
used to determine the correspond:! rig temperature distribution. Fewer 
iteration steps are required, however, if a more refined solution, refer- 
ences 2 and 5, for example, is taken for the first approximation. With 
the initial, values determined by the method of reference 2, four or 
five iterations were required to obtain differences of less than l/2 per- 
cent between the last two iterations (for both the shear stress and the 
heat-transfer rate at the surface) . 

The following procedure was found to be most satisfactory in the 
numerical solution. With the initial values of velocity and temperature 
given for a specified surface temperature and Mach number as a function 
of the independent variable t] and with the variation of viscosity and 
Prandtl number with temperature prescribed, the functions f(q), a(q), 

and J(ti) are determined by numerical integration (the trapezoidal rule 
was found to be very satisfactory) . The velocity distribution is then 
obtained from equation (13) • By makl rig use of the shear function given 
by equation (12) and the previously calculated values of f (tj) , the 
functions P(t]), ^t]), U(tj) , and V(t]) are determined through numerical 

integration and the enthalpy function -8 (t}) is determined from equa- 
tion (17) . The corresponding temperature distribution may be found from 
the enthalpy distribution with the use of enthalpy tables. With these 
new values of velocity and temperature, the steps outlined above are 
repeated and a second approximation to the velocity and temperature 
distribution is determined. This process is continued until the desired 
accuracy is obtained. The procedure for calculating the velocity and 
temperature distribution for the insulated surface is similar to that 
described for the condition of constant surface temperature. The ent halp y 
function in this case, however, is found from equation (15) . 

Experimental values of the thermal properties or values found from 
analytical expressions may be used with equal facility in the computa- 
tions . The equations may be simplified somewhat for the particular case 
in which the Prandtl number is taken as a constant; then, Pr* = 1 and 

the integrating factor p reduces to aP 1 * . For solutions involving a 
constant Prandtl number, however, the method of reference 5 is preferable 
since the solution can be obtained with less calculation. 


BESULTS AND DISCUSSION 


Some friction and heat-transfer characteristics at supersonic Mach 
numbers have been computed by the method described in the prece ding 
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section. Faired experimental values of the viscosity and Prandtl num- 
ber, taken primarily from reference 11, have been used in the computa- 
tions . (The calculations were lim i ted to the temperature range in which 
experimental data were available.) The free-stream temperature was taken 
as 392. 7° R (-67° F), the value at the isothermal region of the MCA 
standard atmosphere, and the corresponding values of the thermal proper- 
ties were taken as : 

c pi = 7*718 Btu/( slug) (deg) 

= 3*058 x 10“7 slugs/(ft) (sec) 
k-L = 3.227 x 10" 6 Btu/ (ft) (sec) (deg) 

Pt-l = 0.73 

The boundary-layer characteristics presented herein include those 
cases which were calculated in reference 6 ; the same thermal data were 
used in both analyses. Some small differences were found between the 
results of reference 6 and those calculated by the method presented 
herein. These differences are due primarily to the poor convergence 
of the solutions in reference 6. Additional calculations have been 
made for those cases by means of the method of this report in order to 
make the results consistent. 

The calculated velocity and temperature profiles for several super- 
sonic stream Mach numbers for the insulated surface are presented in 
figure l(a) , and profiles for T s * = 1, 2, and 1 are presented in fig- 
ures l(b) , 1(c), and l(d) , respectively. (Tabulated values of the 
velocity and temperature distributions and their first derivatives, to 
four decimal places, together with a chart of the thermal properties 
used herein, can be obtained upon request from the National Advisory 
Committee for Aeronautics.) The corresponding friction and heat-transfer 
characteristics are presented in figures 2 to 6 together with the results 
of several other analyses which make use of various approximations for 
the viscosity and Prandtl number . 

Skin friction.- The shear stress at the surface t is given by 


T 



pi u i 2 1— \ 

vs l a i A=o 
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where R is the Reynolds number U^x/v^. Substituting from equation ( 13 ) 
results in the following equation for shear stress: 

pift 2 1 

T = 

{r j(») 

since a(0) = 1. The average skin-friction coefficient c^ is found 
hy integration as 

c^. = — f t dx = — — — (l8) 

i pjlJj 2 * Jo ft J (") 


The skin-friction coefficients found from equation (l8) for the 
insulated surface and for the condition of heat transfer are presented 
in figures 2(a) and 2(h), respectively, together with the results from 
other analyses, for a stream temperature of - 67 ° F. All the analyses 
show substantially the same variation of skin friction with Mach number 
for the insulated surface except for the solution based upon the linear 
viscosity relation |i* = T* for Pr = 1.0; for these conditions the 
skin friction is independent of the Mach number. The present results 
and those of reference 5, which are based upon a constant Prandtl number 
of 0.75 and the use of the Sutherland equation for viscosity, are in 
excellent agreement for both the insulated surface and for the specified 
constant surface temperatures. The differences between the re exalts of 
this paper, those of Moore (ref. 8), and those of Young and Janssen 
(ref. 7) 1 may possibly be attributed in part to some inaccuracies in 
the differential-analyzer calculations, differences in the thermal data, 
and a small difference in stream condition - the calculations from refer- 
ences 7 and 8 are based upon a nominal stream temperature of 400° R. 


-ktn reference 8, experimental values of viscosity are used for tem- 
peratures up to 3,240° R and for Prandtl number, up to 2,000° R; for 
temperatures between 3>240° R and the dissociation temperature the vis- 
cosity was calculated with the use of the power law 

and for temperatures over 2,000° R the Prandtl number was determined from 
calculated values of the thermal properties . Values of the Prandtl num- 
ber from various sources and the Sutherland equation for viscosity were 
employed in the boundary-layer calculations reported in reference 7* 


T 


T=32kO 


v3240 j 
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The agreement of the results based on a power law and a linear viscosity- 
temperature relation with the results of this analysis is not so satis- 
factory as that of the results which make use of the Sutherland equation 
for viscosity. 

It should be noted, however, that the variation of s kin -friction 
coefficient with Mach number is dependent on the stream temperature and 
thus the differences between the results based on a power law and a 
linear viscosity-temperature relation and the results of this analysis 
would depend on the stream temperature at which a comparison was made. 

With a linear viscosity-temperature relation of the form 
H = Constant T*, the skin friction is independent of Mach number; 

taking the constant as unity results in Cf \/r = I.328 for all sur- 
face temperatures and stream Mach numbers irrespective of the value of 
the Prandtl number. The linear viscosity relation = Constant T* 
was utilized in the calculations of reference 2 for a Prandtl number 
of O.72 and the constant was adjusted to give the correct value of the 
viscosity at the surface. This procedure yields a variation of skin 
friction with Mach number for the insulated surface since the recovery 
temperature varies with Mach number; for the condition of constant sur- 
face temperature, the skin friction is independent of the Mach number 
and depends only upon the surface temperature. 


Heat-transfer characteristics .- The heat-transfer rate per unit 
area at the surface q is given by the expression 


l st\ _ _jh. FT m\ 

V a y / y= 0 PrJ V 1 X 2 \FT* 4r l/r|=0 


and from equation (l6) , since 3=1 and 9=0 for q = 0, 


q = - 


^1 1 . / p l^l U l 7 ~ 1 m 2 " ^s 1 j p l^lPl i e " 1 e 


Pr n 


V(co) 


Pih 


v(°°) 


(19) 


where 7 is the ratio of the heat capacities at constant pressure and 
constant volume . 
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As can be seen from equation (19) , tbe heat-transfer rate is pro- 
portional to the difference between an effective enthalpy function -d e 

and the surface enthalpy function d s (or the difference between the 
effective enthalpy i e and the surface enthalpy i s ) . In contrast, 

the heat-transfer rate is approximately proportional to the difference 
between the recovery temperature (or enthalpy) and the surface tempera- 
ture (or enthalpy) when the Prandtl number is taken as a constant. The 
effective enthalpy function and the recovery enthalpy function (the 
value of -Q at the surface for an insulated plate) are presented in 
figure 3 as a function of the stream Mach number. Both -Q r and i9 e 
decrease with increasing Mach number and d e decreases with increasing 
value s of the surface temperature; the functions -3 r and become 

identical for the condition of zero heat transfer. 

The temperature recovery factor r may be determined from -d r 

with the use of enthalpy tables . The recovery factor is presented in 
figure If as a function of Mach number together with the values found 
from other analyses. With the assunption that the Prandtl number is 
independent of temperature, the recovery factor is approximately equal 
to the square root of the Prandtl number for laminar flow for all Mach 
numibers (refs. 2 and 5, for example) irrespective of the viscosity^ 
temperature relation; thus the recovery factor is unity for a Prandtl 
number of 1.0 and the recovery temperature equals the stagnation tempera- 
ture of the stream. When the variation of Prandtl number with temperature 
is taken into account, however, the temperature recovery factor decreases 
substantially with increasing Mach number. The recovery factors of 
reference 7, based upon the use of the Sutherland equation for viscosity 
and experimental values of Prandtl number, are in substantial agreement 
with the results of the present analysis; the recovery factor is approxi- 
mately 0.85 at low Mach numbers and decreases to approximately 0 .76 and 
0.65 at Mach numbers of 5 and ID, respectively. 

Equation ( 19 ) for the heat-transfer rate is considerably more com- 
plicated than those determined from analyses based upon the use of a 
constant Prandtl number. The function i3 e is dependent upon both the 

surface temperature and the stream Mach number for a given stream tempera- 
ture; whereas the recovery factor (a constant for all Mach numbers) takes 
the place of d e when the Prandtl number is taken as a constant . 

The empirical relation 


*e " ^s 


v(°°) (d r - o 


1 35 

= 0.30532 - 0.001289 2.25 - M x - 0.006062 T g * (20) 
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represents the variation of - 8 e with 3 r within 1 percent for Mach 

numbers up to 5 for a stream temperature of - 67 ° F. Values from equa- 
tion ( 20 ) are compared with the calculated values of the quantities in 


figure 5 . Eliminating — — between equations ( 19 ) and (20) leads 

V(“) 


to the following expression for the heat-transfer rate: 


q = -- 7 — ~ . 1 M 2 (0 .30532 - 0.001289 1 2. 25 - mI 1 ^ 5 - 

Pr x |( x 2 

0 .006062 T g *) (* r - $ s ) (21) 

The values of presented in figure 3 and the values of the tempera- 

ture recovery factor presented in figure 4 for Mach numbers between 5 
and 10 were determined from equation (20) with values of -3 e and V(°°) 

calculated for T B * = 2 . With these values of d r the heat-transfer 
rate calculated from equation (21) agrees with those computed from equa- 
tion (19) for T B * = 1, 2, and 4 to within 1 ~ percent for Mach numbers 

from 1 to 10. 

Qomparison of heat-transfer coefficients or Nusselt numbers found 
from different analyses is difficult because of the different tempera- 
ture functions which arise in the expressions for the heat-transfer rate 
A direct comparison of the heat-transfer rates more clearly shows the 
effects of the various approximations on the calculated heat-transfer 
characteristics . The heat-transfer rate determined by the method of 
this analysis is presented in figures 6(a), 6(b), and 6(c) for values 
of T s * = 1, 2, and 4, respectively, together with the results of other 
analyses . 

Figure 6 shows that the local heat-transfer rates calculated with 
the use of the Sutherland equation for viscosity and, a constant Prandtl 
number of 0 .75 are in very good agreement with the results of this 
analysis for all Mach numbers and surface temperatures considered. Like 
wise the results of reference 8, which employ values of the thermal 
properties taken in part from experiment, are in good agreement with 
the results of this analysis for T s * = 1, the only value of surface 

temperature where a comparison could be made . The calculated heat- 
transfer rates which employ more approximate viscosity-temperature 
relations are considerably higher than the values determined by the 
method of this report at very high supersonic Mach numbers; below a 
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Mach number of k or 5, however, the agreement is generally satisfactory 
for a stream temperature of - 67 ° F. 


CONCLUDING REMARKS 


Friction and heat-transfer characteristics of the laminar boundary 
layer on a flat plate have been calculated for Mach numbers from 1 
to 10 for a stream temperature of - 67 ° F with the use of experimental 
values for the heat capacity, conductivity, and viscosity. The conse- 
quent variability of the Prandtl number throughout the boundary layer 
causes the recovery factor to decrease from a value of approximately 
0.85 at low Mach numbers to approximately O .76 and 0.65 at Mach numbers 
of 5 and 10, respectively; whereas the recovery factor is approximately 
equal to the square root of the Prandtl number for all Mach numbers 
■when the Prandtl number is independent of temperature . Moreover, the 
heat-transfer rate is proportional to the difference between an effective 
enthalpy, which is a function of both the surface temperature and stream 
Mach number, and the surface enthalpy when the temperature variation of 
all the thermal, properties is considered; whereas the heat-transfer rate 
is approximately proportional to the difference between the recovery 
enthalpy and surface enthalpy for calculations based upon a constant 
Prandtl number. Because of the excellent agreement between the results 
of the present analysis and calculations based upon the use of the 
Sutherland equation for viscosity and a Prandtl number of 0.75, it 
appears that the exact variation of the thermal properties of the fluid 
need not be considered in the calculation of shin friction or heat 
transfer. More approximate viscosity-temperature relations and the use 
of a Prandtl number of unity introduce significant errors in the calcu- 
lated friction and heat-transfer characteristics at Mach numbers above It- 
or 5 for a stream temperature of - 67 ° F. 


• Langley Aeronautical Laboratory, 

National Advisory Co mmi ttee for Aeronautics, 
Langley Field, Va., January 23, 1953* 
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Figure 1.- Velocity and temperature profiles. 
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Figure 1.- Continued. 
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“r 

i r Extrapolated!, from 
equation (zOj for T g *=2 


Free- stream Mach number, 


Present analyile 

Sutherland eq. ; Pr=0.75 (ref. 5) 

li'-oT*; Pr-O .72 (ref. 2) 

Extrapolated, from eq.^c} for T g *=2 

xV Sutherland eq., Pr variable (ref. 7) 



Free-stream Haoh number, 


Figure 4.- Variation of temperature recovery factor with, stream Mach n umb er 
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Figure 6.- Variation of local heat-transfer rate with stream Mach number. 
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(b) T * = 2. 
v ' s 

Figure 6.- Continued 
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Free-stream Maoh number, M 


(c) T * = k. 
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Figure 6 . - Concluded . 


NACA-Langley - 8-31-64 - UO 



